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The phase transition occurring in a square 2-D spin lattice governed by an anisotropic Heisenberg 
Hamiltonian has been studied according to two recently proposed methods. The first one, the 
Dressed Cluster Method, provides excellent evaluations of the cohesive energy, the discontinuity of 
its derivative around the critical (isotropic) value of the anisotropy parameter confirms the first- 
order character of the phase transition. Nevertheless the method introduces two distinct reference 
functions (either Neel or XY) which may in principle force the discontinuity. The Real Space 
Renormalization Group with Effective Interactions does not reach the same numerical accuracy but 
it does not introduce a reference function and the phase transition appears qualitatively as due to 
the existence of two domains, with specific fixed points. The method confirms the dependence of 
the spin gap on the anisotropy parameter occurring in the Heisenberg-Ising domain. 



I. INTRODUCTION 

The study of spin or electron lattices, even when they 
are governed by simple model Hamiltonians, requires in 
general approximate methods in order to obtain reliable 
estimates of the cohesive energy, of the excitation gap, 
of the spatial correlation, etc... The treatment of phase 
transitions is a special challenge for approximate methods 
since it is in general not easy to identify the values of the 
interactions at the critical points, the nature of the phase 
transition, as well as the behavior of the properties on 
both sides of the phase transition. The purpose of the 
present work is to compare the abilities of two methods 
recently developed by the authors to study a first-order 
phase transition. 

Despite its rather formal character the spin ^ anisotropic 
Heisenberg Hamiltonian on an infinite 2-D square lattice 
may be used as an excellent model problem to test the 
ability of a theoretical method to treat a phase-transition 
phenomenon. This Hamiltonian is given by 

H = JJ2 (SfSJ + SfS] + \StS^), (1) 

(ij) 



where a and (3 are the usual spin up and spin down func- 
tions, and those of the other sublattice bear 

X^{a~- /3)/\/2. (3) 

If one works in the basis of (X,Y) functions instead 
of {a, (3) ones, this XY polarized function will appear 
as the leading configuration for the — 1 < A < 1 do- 
main. At A = 1 (isotropic Hamiltonian) a transition 
to an Ising-like phase occurs. Actually for A ^ oo 
the Hamiltonian becomes an Ising Hamiltonian and the 
ground state becomes the Neel fully spin-alternate func- 
tion $0 = a(3a(3..., which is also the leading configura- 
tion for A > 1. Early Quantum Monte Carlo (QMC) 
calculations^ suggested, although with some imprecision, 
that this transition is of first-order type. More recent 
and more accurate calculations (see for instancrfS) have 
confirmed its first-order character. One may also quote 
elaborate Coupled Cluster (CC) calculationsSAifl which 
start from either a planar-like function or the Neel wave 
function as reference function $0 and assume an expo- 
nential form of the wave operator 

I*) = expS\<^o), (4) 



where runs over all pairs of nearest neighbor 

sites. This 2-D square lattice model has no exact 
solution and has therefore been the subject of numerous 

Calculation9ii2i2iiiSiLLLSiifliiLiSii2iiiiiSiiLiLiSii2i2SiSiiS2i2S 

in the recent past, which employ either ana- 
lytic expansionSfi^'i^'iS'i^ or numerical algorithms 
such as Coupled Cluster approaches jSiSiiSiii exact 
diagonalizationsiiii^ and Quantum Monte Carlo 
calculationsii .22.23,24,2.5,2^,27 At A = -1 a first- 
order transition takes place between the ferromagnetic 
phase and a planar-like phase in which the spins in the 
ground state wave function lie in the XY plane. This 
so-called XY polarized function is such that the sites of 
one sublattice bear 



where S is restricted to a certain number of local many- 
body operators (up to 6-body operators). The results 
agree very well with those of QMC calculations in the two 
24.2ieg^.§ps around A = 1, each region being treated using the 



(2) 



relevant reference. Although the authors do not conclude 
expHcitly, the results support the first-order character of 
the phase transition at A = 1. The extent of the domain 
of bi-stability is more difficult to assess since it seems to 
depend on the sophistication of the wave-operator. The 
present work studies the same problem using two new 
methods have different characteristics. The methods em- 
ployed hereafter 

- the Dressed Cluster Method (DCM?!) uses, as do 
the Coupled Cluster expansion (CC) and pertur- 
bative approaches, a single reference wave function 
$0, which will be either the Neel function or the XY 
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polarized configuration. In DCM this wave func- 
tion is used as a bath in which a finite cluster is em- 
bedded and treated exactly. Then the configuration 
interaction matrix relative to the cluster is dressed 
under the effect of excitations occurring around 
the cluster, the amplitudes of which are transfered 
from the amplitudes of similar excitations within 
the cluster. This approach will be shown to give 
extremely accurate results, very close to the best 
Quantum Monte Carlo calculations of the cohesive 
energy and confirm the first-order character of the 
phase transition but, as well as the CC method, it 
suffers from the prejudice introduced by the discon- 
tinuity of the reference function <i>o . 

- the Real Space Renormalization Group with Ef- 
fective Interactions (RSRG-El22i) is an improve- 
ment of the RSRG method originally proposed by 
Wilson^i It proceeds through the same reduction 
of the Hilbert space by considering fragments (or 
blocks) of the lattice, and a reduction of the Fock 
space for these blocks to a few states of lowest 
energy. But it extracts effective interactions be- 
tween the blocks through the exact diagonaliza- 
tion of dimers of blocks. The knowledge of the 
exact spectrum of the dimers enables one to de- 
fine, using the theory of effective Hamiltonians pro- 
posed by Bloch(21 inter-block effective interactions. 
The method is iterative, it is repeated to blocks 
of blocks, etc... until it reaches fixed points of the 
problem. The method provides at a very low cost 
reasonable estimates of the cohesive energy of 1- 
D or 2-D spin lattices. It does not introduce any 
reference function, it is therefore in principle con- 
tinuous on both sides of the critical value of the 
parameter. However the method leads to two dis- 
tinct fixed points for the A < 1 and A > 1 domains. 
The iterations result in a discontinuity of the cohe- 
sive energy derivative. The method also shows the 
appearance of an excitation gap for A > 1. 



II. DRESSED CLUSTER METHOD 

Let us summarize the main points of the Dressed Clus- 
ter Method : 

- one first defines a single-determinantal reference 
function $0 on the infinite lattice, namely the Neel 
or the XY function. For sake of simplicity, the 
method will be presented here using only the Neel 
function in the a(3 representation 

$o = n2*(2« + l), (5) 

i 

- one considers a 2-D square finite cluster of N sites 
which divides the atoms in two subsets, internal 



and external, so that the reference function can be 
written as 

$0 = (6) 

- the model space S is spanned by the determinants 
obtained from $0 by all possible excitation pro- 
cesses which only concern atoms within the 
cluster 

s = {$,} = {$g-*.r+$™*}. (7) 

Let Ps be the projector onto this model space. The di- 
mension of the full Configuration Interaction (CI) space 
is equal to that of the isolated cluster. Nevertheless the 
diagonal elements of the matrix PgHPs differ from those 
of the isolated cluster CI matrix under the effect of the 
embedding, i.e., the energy of each determinant is shifted 
by a quantity J; per alternating bond 1, at the frontier, 
the determinants $i in the lattice problem interact only 
with the outer-space determinants D^'^i obtained from 
<i>i by a spin exchange Df on the external bond /. Re- 
placing for simplicity the determinants by their index 
i, the eigenequation for fine i is 

j&s.j^i I ext 

The last summation must be evaluated through a proper 
estimate of the coefficients Cn+. . These coefficients are 
approximated to the product of the coefficients of the 
determinants <I>i by environment-dependent amplitudes 
di^i characteristic of the excitations Df on $i. 

= cal,^. (9) 

These amplitudes are extracted from the knowledge of 
the CI wave function of the embedded cluster. 
In order to be more explicit, let us consider a determinant 
<i>i (cf Fig.l). The cluster is delimited by a continuous- 
line box and is embedded in the Neel function. Bonds 
involved in the excitations from $0 to <I>i appear with 
thick lines. The elementary excitation D'^ on an external 
bond 1 (indicated by a dashed line) leads to a determi- 
nant which interacts with <I>i through an exchange 
integral J;. The excitation amplitudes di^i depend on the 
environment of the bond 1 (the largest considered envi- 
ronment is indicated by a dashed-line box) and are taken 
as 

dl,r = (10) 

where r is a translation from the external bond I to the 
outermost equivalent bond ^ — r of the cluster (which is 
indicated by a continuous line) and S 5 is such that 
the environment of the bond Z — r in <I>j has the maxi- 
mum resemblance with the environment of bond Z in $j . 
One must notice that, in some cases, it is necessary to 
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FIG. 1: Dressed Cluster Method: schematic view of the prin- 
ciple of the Configuration Interaction dressing under the effect 
of the spin exchanges around the cluster, pictured by a full 
line box. The upper part identifies the cluster and a determi- 
nant embedded in a Neel environment, as well as the outer 
bond 1 on which a spin exchange will be performed. The lower 
part pictures the two determinants from which the amplitude 
(Eq. 10) will be extracted. 



restore the right spin = of the translated determi- 
nants by changing the spins of the atoms furthest from 
the bond 1, in order to obtain the most relevant infor- 
mation from the CI wave function. Finally the quantity 

E; ext H.^D+ fio+i can be replaced by 

( Jldl,^)C^. (11) 

I ext 

This summation can be delt with as a diagonal energy 
shift (dressing) 

A,, - ( ^ Jid,,,;), (12) 

/ ext 
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FIG. 2: Cohesive energy as a function of the anisotropy pa- 
rameter A. (+) DCM (from a 16-site cluster) with the Neel 
reference function, (□) DCM (from a 16-site cluster) with the 
XY polarized reference function. 

and the corresponding dressing operator A 

A = ^|ci>,)A,,(<i>,|, (13) 

Eq. 9 insures the translational invariance; if the deter- 
minant D'^^i is identical through a translation T to one 
of the determinants belonging to S i.e., if D^^i = T 
then C^+$. = C**;,- This estimation of C^+j leads to 
an important simplification : the effect of excitations on 
bonds I which are far from the fragment (by more than 
the cluster size) is approximated to be identical for all 
determinants <i>i and only shifts the diagonal elements of 
the CI matrix by the same amount. It has consequently 
no effect on the eigenvectors of the dressed CI matrix 
Ps{H + A)Ps and can be omitted. Since the dressing de- 
pends on the eigenvector the procedure must be repeated 
to self-consistency. One may say that the DC method im- 
plies many-body operators, up to the number of atoms in 
the cluster. It does not proceed to a strict exponential- 
ization of the wave function but it employs ratios of co- 
efficients to transfer information from the internal CI to 
take into account the effect of elementary excitations on 
the external bonds. Through the environmental depen- 
dence of these elementary excitation amplitudes, many- 
body effects are introduced. The relation with a Coupled 
Cluster expansion of the wave function2Si52i2ii2i has been 
discussed in ref. 13. The accuracy of the DC method has 
been illustrated on 1-D electron and spin (frustrated and 
non-frustrated) lattices. It has also been applied to the 
study of the lowest excitation energies as functions of the 
bond alternation in the 1-D spin chainiS^: The DC method 
is now applied to the 2-D square spin lattice using a 4 x 4 
cluster and starting from both the Neel function and the 
XY polarized function as reference <i>o- The computed 
cohesive energy as a function of A is pictured in Fig. 
2, where the two branches, obtained from the XY and 
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FIG. 3: Cohesive energy around the isotropic point, (o) 
QMCf21 (□) DCM, (+) RSRG 



Neel functions respectively, appear clearly as crossing in 
A = 1. One observes the existence of a continuation of 
the Neel-generated solution in the 0.4 < A < 1 and of 
the XY-generated solution in the 1 < A < 1.5. This may 
be seen as the indication of metastable states around the 
critical A = 1 value, as expected for a first-order phase 
transition. 

The quality of the DCM results has to be assessed by 
comparison with accurate analytical or numerical calcu- 
lations. Table I and Fig. 3 report such comparisons. 
For A = 1 our estimate — 0.66928J coincides to 10~^ 
with the most accurate QM02Si2£i21 value -0.66944J. It 
may be interesting to compare with CCM result o^i^i^° 
which are — 0.6670J when introducing 6-body opera- 
tors, and — 0.66817J when introducing 8-body operators. 
The difference indicates the importance of many-body 
operators, and the slow convergence in this expansion. 
The 3'''^-order spin-wave gives —0.6700 J and a plaquette 
expansioiJi — 0.6691J. 

The agreement of our DCM values with QMC calcu- 
lations is similar for A ^ 1. For A = we obtain 
-0.5489 J, similar to the result of Lin et al2i -0.54882 J, 
or for A = 0.6 (DCM = -0.61094J, QMC = -0.60958J). 
Fig. 3 shows the near identity of our results with those 
of Lin et al in the whole < A < 1 domain. For 
A > 1, the agreement in similar, as may be seen from 
Fig. 3 and Table I. For instance we obtain nearly iden- 
tical values for A = 2 (DCM = -1.08329 J, QMC = 
— 1.08220J). Notice that we have no convergence prob- 
lem when A ^ l"*", while cohesive energies could not be 
obtain in the 1 < A < 1.09 domain in ref. 27. It is 
clear that DCM represents, in view of its low cost, an 
interesting alternative to QMC. 



TABLE I: Cohesive energy of the anisotropice 2-D lattice. 



_A DCM RSRG QMO^^ 

-0.54890 -0.53966 -0.54882 
0.6 -0.61094 -0.60260 -0.60958 

1 -0.66972 -0.66615 -0.66944 
1.2 -0.73961 -0.73072 -0.73920 

2 -1.08329 -1.07849 -1.08220 



III. REAL SPACE RENORMALIZATION 
GROUP WITH EFFECTIVE INTERACTIONS 

A. Method 

The Real Space Renormalization Group proposed by 
Wilson essentially consists in an iterative truncation of 
the Hilbert space. The method proceeds through the def- 
inition of blocks of N sites, periodisable fragments of the 
periodic lattices and the research of the (lowest) eigen- 
states of the Hamiltonian relative to these blocks. For a 
block I, and the corresponding Hamiltonian Hj 

Hi(l)Kj ^ Sk,i4>kj- (14) 

One shall retain a few (let say m) eigenstates of Hi . Then 
one will consider a block of blocks (1...I...J...N), and one 
will approach the wave function for this superblock by 
working in a truncated Hilbert space constituted of all 
products of the m eigenstates kept for each block. 

W ^K,i, K = l,m. (15) 

I=l.N 

Then the process can be repeated, till convergence. If the 
blocks and the sets of selected eigenstates are properly 
defined the problem at each iteration may keep its for- 
mal structure, while the interactions between the super- 
super... sites change along the iterations. One then 
reaches in a certain number of steps a fixed point of the 
problem. 

This idea is extremely elegant. However the attempts to 
use it as a practical numerical tool for the study of peri- 
odic lattices (of either spins or electrons) were extremely 
discouraging. And the method was abandoned, although 
it gave birth to a deeply different formalism, namely 
the Density Matrix Renormalization Group, which is 
extremely performant, but limited to the treatment of 
(quasi) 1-D systems. 

The failure of the RSRG method is due to the simple 
truncation of the Hilbert space and the total neglect 
of the non-selected eigenstates of the blocks. Rather 
than trying to treat the effect of the non selected states 
in a 2"''-order perturbative approachjSI two of the au- 
thors have suggested to define effective-interactions be- 
tween adjacent blocks A and B by solving exactly the 
Schrodinger equation for the AB dimer, and by mak- 
ing use of the Bloch's theory of effective Hamiltonians. 
We shall not repeat here the formalism, given in ref. 
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29, which leads to a modified RSRG formalism, called 
RSRG-EI (EI = Effective Interactions). The first test 
applications of the method were quite encouraging. We 
simply make explicit hereafter the specification of the 
method in its simplest version for the study of a square 
spin lattice. 

The method consists in considering a square (3 x 3) block 
of 9 atoms. Its ground state is a doublet with Sz = ±1/2 
and it is the only state kept hereafter. Let call a and o 
the Sz = 1/2 and = —1/2 degenerate doublet ground 
states of the block A. The block can therefore be seen as a 
super-spin. In order to establish the effective interactions 
between the ground states of adjacent blocks, one treats 
exactly the 18 (3 x 6)-site superblock AB. One wants 
to establish the effective energies of and interactions be- 
tween the four products of ground state wave functions 
which define a model space ab, ah, ah, ah. Diagonal- 
izing the exact Hamiltonian for the AB superblock one 
may identify the eigenstates ^'^(S'z = 1), ^^{S^ — — 1), 
'i!^{Sz = 0), ^'^(S'z = 0) which have the largest projec- 
tions on the model space, and their energies E'j'+ = E^^ , 
Erpo and Ego. The three energies can be seen as the 
eigenvalues of a new anisotropic Hamiltonian 

+Ea + Eb + AEab- (16) 

Hence 

Et+=Ea + Eb + AEab, (17) 



Eto = -- J(i)a(i) + - J(i) +Ea + Eb + /^ab, (18) 



Ego = -- J(i)a(i) - - J(i) +Ea + Eb + AEab. (19) 
From which one obtains 

J^^^ ^ Ego - Eto, (20) 



j(i)a(i) = 2Et+ - Eto - Ego. (21) 

These equations define a new anisotropic Heisenberg 
Hamiltonian between blocks. The process may be re- 
peated, treating a block of 9 blocks and a superblock of 
18 blocks, till convergence is achieved. 

B. Results 

The qualitative key points in that problem are the facts 
that 




FIG. 4: RSRG study: anisotropic parameter "I"'^^ after the 
first iteration, as a function of the initial anisotropic param- 
eter $ (Eq. 20). The stairs illustrate the convergence of the 
iterative procedure to the fixed points, Ising on the right side, 
XY on the left side 

(i) for A = 1, A'^-' — 1, the problem remains isotropic 

(ii) for A > 1, A*^^-' > A, the anisotropy is increased in 
the direction of an Ising problem 

(iii) for A < 1, A'^^-' < A, the anisotropy increases in the 
opposite direction towards a pure XY problem. 

For graphical purposes the anisotropic Hamiltonian may 
been written as 

H = J[[SzSz)sm $ + (5,5, + SySy) COS $]. (22) 

The isotropic case corresponds to $ = — 7r/4, the XY 
problem to <i> = 7r/2, the Ising situation to $ = — 7r/2. 
On sees that A = tan<I>. Fig. 4 reports the evolution 
of $(1) as a function of $. The iterative process, start- 
ing from as new value <i> leads to a new anisotropy angle 
<i>i = <i>(i)($). The second step leads to <i>2 = <i>(^)(<i>i), 
etc... The qualitative nature of the phase transition ap- 
pears dramatically. Starting from $ > 7r/4, in- 
creases rapidly. As seen from Fig. 3 the process con- 
verges in a few steps to the = 7r/2 fixed point, i.e, 
to an Ising problem. Oppositely, starting from $ < 7r/4, 
<i>i decreases. The fixed point on that side A < 1 is 
the pure XY problem (A = -1, $ = -7r/4). But the 
curve — /($) is tangent to the fine of slope one 
$(1) — $ for $ — — 7r/2. Hence the fixed point is in prin- 
ciple reached in an infinite number of steps. 
The quantity J'-^-' is significantly lower than one for A < 1 
and tends to zero when A tends to -1. It increases with A 
but remains finite in the region A > 1. Fig. 3 reports the 
RSRG-EI calculated cohesive energy. For A = 1, as al- 
ready reported elsewhere^ the RSRG-EI cohesive energy 
is Ecoh = — 0.666155J. This value is in slightly poorer 
agreement with the best QMC value — 0.66934J than the 
previously reported DCM value, but it is obtained at a 
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0.1 - 




FIG. 5: Appearance of the gap in the A > 1 phase, as cal- 
culated from the RSRG-EI. The full line is proportional to 

(A^ - 1)1/^ 



much lower cost. The underestimation of the cohesive en- 
ergy by the RSRG-EI method is systematic but it never 
exceeds 2% (cf Table I and Fig. 3). We have carefully 
checked the existence of a discontinuity of the slopes of 
the curve Ecoh = /(A) around A = 1. This discontinu- 
ity clearly appears from the insert of Fig. 3. The slope 
{dE/dX)^^^+ between A = 1.02459 and A = 1 is 0.32 in 
QMC and 0.26 in RSRG, on the A < 1 side the slope from 
QMC is 0.175 (between A = 0.97 and A = 1), which the 
slope from RSRG is 0.20 (between A = 0.95 and A = 1), 
0.21 (between A = 0.99 and A = 1). Although weaker 
than the interpolated estimates from QMC, the disconti- 
nuity of the slope predicted from RSRG-EI is clear. The 
existence of a discontinuity was not a priori evident since 
the quantities J^^^ and A^^^ are continuous functions of A. 
The discontinuity comes from the fact that the iterations 
tend to different fixed points for A > 1 and A < 1. 
Actually the method is also able to explain the absence 
of a gap for A < 1 and of the existence of a gap for A > 1. 
For A < 1, since one must repeat an infinite number of 
iterations with decreasing values of J*-^-* , the lowest states 
are degenerate. In the A > 1 domain, the system will be 
gapped since the process converges in a finite number of 
steps, with finite values of J. Fig. 5 reports the calcu- 
lated gap for A slightly larger than 1. We have checked 
the behaviors of the gap as a function of A. Spin- wave 
theory predicts that it should follow the law 

AS = 2(A2- 1)1/2. (23) 

Previous numerical works^i have shown that the excita- 
tion energies are significantly lower, by a factor close to 
0.5. Fig. 5 have used an interpolation 0.86634(A2 - l)i/2 
which fits well our calculated values. 



IV. CONCLUSION 

The present paper studies the behavior of a 2-D square 
spin lattice obeying an anisotropic Heisenberg Hamilto- 
nian. Since it presents a phase transition, this problem 
can be seen as a convenient test to compare the abilities 
of the methods available for the treatment of 2-D (or 3-D) 
lattices. One may subdivide the methods in two groups. 

- methods which rely on (or require the introduction 
of) a simple zero order wave function. This wave 
function may be perturbed, or considered as the ref- 
erence function for a Coupled Cluster expansion (i.e 
an exponential development of the wave operator) . 
In such a case different zero-order or reference wave 
function will be used for the two different phases. 
This choice of two distinct references may be seen 
as forcing the phase transition and presents the risk 
to impose artefactual discontinuities. The here em- 
ployed Dressed Cluster Method only uses the refer- 
ence function as a bath around a finite cluster, but 
it is subject to the same criticism. 

- prejudiceless methods which do not bias the treat- 
ment by introducing reference wave functions. 
Among them one may quote finite cluster exact 
diagonalization, followed by extrapolations on the 
cluster size. For 2-D systems extrapolations are 
quite difficult to perform. Quantum Monte Carlo 
calculations require both statistics and extrapola- 
tion and the error bars may prevent a clear assess- 
ment concerning the nature of the phase transition, 
when for instance the change of the slope of the co- 
hesive energy as a function of the internal parame- 
ter is small. Recent progresses have reduced these 
uncertainties. 

The excellent agreement of the DCM results with the best 
QMC calculations for A = 1 gives confidence in the accu- 
racy of the calculated dependence of the cohesive energy 
on the anisotropy parameter and assesses the first-order 
character of the phase transition. 

The RSRG-EI treatment does not enable one to reach 
such a numerical accuracy but it presents several advan- 
tages 

- it does not introduce the bias of a reference function 

- it visualizes qualitatively the phase transition in 
terms of a critical value of the parameter separating 
two domains with their specific fixed points 

- it offers a simple understanding of the gapless- 
gapped character of the two phases. 

The philosophy of the RSRG method is responsible for 
this qualitative and pictorial advantage. The introduc- 
tion of effective interactions adds a numerical improve- 
ment to this conceptual tool. Of course, as shown for 1- 
D lattices, the quantitative performance of the RSRG-EI 
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treatment is much better when it is possible to extrap- 
olate its results with respect to the size of the blocks. 
This is not possible for the present time for 2-D lattices, 
since the next size of a square block would be 25 (which 
would require the exact treatment of a 50-site problem 



for the superblock) . But the accuracy of the results from 
9-site blocks is surprisingly good and the elegance of the 
method suggests to consider it as an excellent exploratory 
tool. 
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